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Characterization of SINR Region for Multiple 
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Systems 
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Abstract 

This paper considers a wireless communication network consisting of multiple interfering multicast 
sessions. Different from a unicast system where each transmitter has only one receiver, in a multicast 
system, each transmitter has multiple receivers. It is a well known result for wireless unicast systems 
that the feasibility of an signal-to-interference-plus-noise power ratio (SINR) without power constraint 
is decided by the Perron-Frobenius eigenvalue of a nonnegative matrix. We generalize this result and 
propose necessary and sufficient conditions for the feasibility of an SINR in a wireless multicast system 
with and without power constraint. The feasible SINR region as well as its geometric properties are 
studied. Besides, an iterative algorithm is proposed which can efficiently check the feasibility condition 
and compute the boundary points of the feasible SINR region. 

Index Terms 

Wireless multicast system, power control, signal-to-interference-plus-noise power ratio (SINR), 

SINR feasibility, SINR region, Perron-Frobenius eigenvalue. 

I. Introduction 

In wireless communication systems, interference is an inherent phenomenon. Due to the 
broadcast nature of wireless channels, interference arises whenever multiple transmitter-receiver 
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pairs are active concurrently in the same frequency band, and each receiver is only interested 
in retrieving information from its own transmitter. For a particular receiver, the received signal 
is a superposition of its desired signal, interfering signals and background noise. SINR, defined 
as the power of desired signal divided by the sum of the power of interfering signals and the 
power of noise, is a widely used performance measure for wireless communication systems. It is 
analogous to signal-to-noise ratio (SNR) used for single user communication, which has clearly 
understood implication on the bit error rate (BER) and capacity for additive white Gaussian 
noise (AWGN) channels. Using SINR as a surrogate for BER and capacity implicitly assumes 
that the interference is an AWGN. Although there are limitations of this assumption, as reported 
in 0.0. the importance of SINR has never been doubted. 

For a system consisting of multiple point-to-point communication sessions, also referred to as 
unicast system, the SINRs of all receivers form a vector. The feasible SINR region includes all 
the SINR vectors that can be achieved by some transmission powers. The geometric properties 
of feasible SINR region has been studied in (3]-[5j. Reference [3) proves that in the case of 
unlimited transmission power, the feasible SINR region is log-convex. In [4j, it is shown that 
under a total power constraint, the infeasible SINR region is not convex. Reference [5| considers 
a system with only three transmitter-receiver pairs without power constraint, and shows that the 
feasible SINR region is concave. It also provides certain technical conditions under which a 
concavity result for systems with a general number of users is established. In [[6], for the cases 
that the transmission powers are subject to arbitrary linear constraints, a mathematical expression 
for the boundary points of the SINR region is obtained. 

In this paper, we consider the feasible SINR region for systems consisting of multiple point- 
to-multipoint communication sessions, also referred to as multicast system. Multicast enables 
data to be delivered from a source node to multiple destination nodes. Practical examples of 
such configurations include cellular networks and two-way relay networks. In cellular networks, 
a base station multicasts a file to multiple mobile devices that request the file at the same time 
(7). In two-way relay networks, when network coding is applied, a relay multicasts the coded 
packets to two sink nodes [8j. The power control and scheduling for wireless multicast systems 
have been studied in [j9j|—[ 1_2] and the references therein. All these works aim to either minimize 
the system power or maximize the system throughput, subject to the constraint that the SINR 
of all receivers are larger than a given threshold. The feasibilities of the problems, however, are 



unknown. 


For a wireless unicast system, the feasibility of an SINR vector without power constraint 


is determined by the Perron-Frobenius eigenvalue of a nonnegative matrix [13], [14]. In this 
paper, we generalize this result to a wireless multicast system. We first propose a necessary and 
sufficient condition under which an SINR is achievable without power limitation. Based on this 
condition, we figure out the feasible SINR region by giving its boundary points. The approach 
is to find the farthest point of the feasible SINR region from the origin in a given direction. 
An iterative algorithm is proposed to find the farthest point, which is also a distributed power 


control algorithm to solve the power balancing problem [13] aimed to maximize the minimal 
SINR of all receivers. 

Then we analyse the geometric properties of the feasible and infeasible SINR regions. It is 
found that the feasible SINR region of a multicast system is in fact the intersection of the feasible 
SINR regions of all its embedded unicast systems. Based on the results in [3|-[|5] for unicast 
systems, we show that the feasible SINR region of a multicast system is log-convex, and the 
infeasible SINR region of a multicast system with two multicast sessions is convex. We also 
show by an example that, the convexity property of the infeasible SINR region does not hold 
for a general multicast system with more than two multicast sessions. 

Later, the necessary and sufficient condition for the feasibility of an SINR in a multicast system 
is extended to include linear constraints on the power. This result generalizes the results in [6| 
for unicast systems. Besides, in [6j, the zero-outage SINR region for a time-varying system is 
also considered, where the channel gains are selected from a finite set. Suppose the transmission 
powers are not allowed to vary with the channel gains, we establish a reduction that maps any 
instance of the zero-outage SINR problem in a time-varying unicast system to a corresponding 
instance of the feasible SINR problem in a time-invariant multicast system. The idea is to regard 
the multiple receivers in a multicast session as an identical receiver that can experience a finite 
set of channel gains. 

The rest of the paper is organized as follows: In Section |ITJ the system model and problem 
formulation are presented. The necessary and sufficient condition on the feasibility of an SINR 
vector is provided in Section 1111} Section ITV] gives the characterizations of the SINR region and 


proposes an iterative algorithm. Section [V] studies the geometric properties of the feasible SINR 
region. Section [VT] extends the study to include power constraints. Finally, the paper is concluded 







Fig. 1: Example of a multicast network consisting of two multicast sessions. Transmitter R 
wants to transmit data to both R\ and R\. Transmitter T 2 wants to transmit data to both R\ and 
f? 2 - Their transmitted signals interfere with each other. The solid lines represent intended links 
and the dashed lines represent interfering links. 


in Section IVII1 

Notation : The following notations are used throughout this paper. Vectors are denoted in bold 
small letter, e.g., x, with their /1h entry denoted by x % . They are regarded as column vectors unless 
stated otherwise. Matrices are denoted by bold capitalized letters, e.g., X, with X l3 denoting 
the (i, j)th entry. Vector and matrix inequalities are component-wise inequalities, e.g., x > y if 
Xj > r/i for all z; X > Y if X Vj > Y r] for all i and j. The cardinality of a set is denoted by 
“| • |”. The Euclidean norm of a vector is denote by “|| • ||” . The transpose of a vector or matrix 
is denoted by (-) r . I represents an identity matrix with compatible size. 0 represents a vector 
with compatible size whose entries are all zero. 

II. System Model 

Consider a general wireless communication network consisting of N multicast sessions. The 
N transmitters are denoted by R for i = 1,..., N. Each R wants to multicast common data 
packets to K, receivers, denoted by R** for k t = 1,..., K,. Without loss of generality assume 
K, > 1. If K, = 1 for all i, the scenario reduces to the unicast case. The total number of 
receivers in the system is K = R%- Define K, j = {1, 2,..., KR for i = 1 ,N. Fig. [I] 

illustrates an example of such network. Let p, be the transmission power of transmitter R and 
p = \pi,... ,pn] T ■ The channel gain between T, and /?,■'’ is denoted by g ^ . All the multicast 
sessions share the same channel and thus interfere with each other. We assume that interference 


caused by simultaneous transmissions is treated as AWGN with variance identical to the received 






power. The SINR of receiver R,f‘ is given by 


7*‘(pj = 


g r k pu Pi 


( 1 ) 


Y, j #9 r *it.P3 + v 2 ' 

where a 2 is the variance of the background noise and without loss of generality, it is assumed 
to be identical for all receivers. We define the SINR of the i-th multicast session as 


7i( P ) = min (t?(p)} ■ 

The SINR vector of the system is 

r (p) = [ 7 i( p)5 T2(p) 5 ■ ■ ■ j 7 tv(p)]- 

In this paper, we analyze the feasible SINR region of a multicast system, that is, 

T = (r(p) e R n : P > o, P e R N }. 

Proposition 1. Given a vector p 6 R ,v . There exists a power vector p* > 0 such that T(p*) = p, 
if and only if there exits a power vector p' > 0 such that T(p' ) > V- 


Proof: The “only if” part is trivial and we show the “if” part. Suppose 7 ;(p') > //, for 
some i. Fix such an 1 . Since ydpj is monotonically decreasing as p, is decreasing, we can find 
a 0 < pf ] < p' and let p (1 > = [pi,... ■ ■ ■ iPnY such that 7i(p (1) ) = On 

the other hand, since 7 j(p) for j f 1 is monotonically increasing as p, is decreasing, we have 
7 ? (p : l ) ) > pj. By keeping on decreasing the power of transmitters that achieve higer SINR than 
p, we obtain a sequence pf\pf \ ■ ■ ■ ,pf \ ■ • • for each i — 1,..., N. It can be seen that these 
sequences are monotonically decreasing and lower bounded by zero, so they are convergent. 
Denote the limit point by p*. For any arbitrarily small 5 > 0 and for all i, since 7,(p) is 
continuous with p, there exists a sufficiently large T, when t > T, | 7 ,;(p*) - 7 ?:(p (t) )| < 6. 
Meanwhile, since 7i(p < - t, - ) ) = Pi for some f > T, we have |7i(p*) — Pi\ < 5. Therefore 7i(p*) = 
Pi for alii. ■ 

By Proposition [I] we say that an SINR vector p, = [pi,..., p N \ is feasible if and only if there 
exists a power vector p > 0 such that 
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In matrix form, it is 


A(ai)p > n(/i), 
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where 
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Each row of A(/x) corresponds to a receiver. For the convenience of discussion, we use af' e 


pTV 


to denote the row of A(/z) that corresponds to receiver i?**. As in the form Q, the feasibility 


of fx can be checked through linear programming [151. However, in a different way, we propose 
a necessary and sufficient condition on the feasibility, which generalizes the Perron-Frobenius 
eigenvalue criteria for unicast systems (square matrices). This condition is used to explicitly 
characterize the feasible SINR region T and to prove some geometric properties of it. Before 
further discussion, we given some definitions. 

Define set 
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: ki G /Q for i — 1,..., N >. 






























Notice that, for each G £ G{pt), only one receiver is involved for each transmitter, which is a 
unicast scenario. So G(fJt) is the set including all the embedded unicast systems and its size is 
Ki. Considering the example in Fig. [l| the four embedded unicast systems are 



In subsequent discussion, we also use k = (ki, k 2 , ..., k N ) to specify a G £ G(n). Let n G = 
[n^ 1 , tt- 2 2 , • • •, n'^ f ] T denote the noise vector with entries of n(/z) that correspond to the receivers 
in G. For the simplicity of notation, we sometimes drop the argument /i of A, n and Q when 
the context is clear. 


Definition 1. [16] A matrix X is called nonnegative if X > 0. A nonnegative square matrix X 
is irreducible if for every pair (■ i,j ) of its index set, there exists a positive integer n = n(i,j ) 
such that X-'- ) > 0 , where X i j " } is the ( i,j)th entry of X''. 


Definition 2. [16] Let X he an irreducible nonnegative square matrix. The Perron-Frobenius 
eigenvalue of X is the maximum of the absolute value of eigenvalues of X, and is denoted by 
A(X). 

Let 1 = [1,..., 1] be a vector whose components are all 1. For each G £ Q{ 1), I — G is the 
normalized interference link gain matrix of the corresponding embedded unicast system. 


Definition 3. A multicast system is called irreducible if and only if the matrices I — G for 
G £ (/(l) are all irreducible. 


It needs to be mentioned that if a multicast system is irreducible, as long as p > 0, the 
matrices I — G for G £ G(fi) are all irreducible. Throughout the paper, we assume that the 
multicast system is irreducible. 


III. Feasibility Condition for SINR 

In this section, we give a necessary and sufficient condition for the feasibility of an SINR 
vector in a wireless multicast system. Recall that in a wireless unicast system, the following 
theorem from [13], [17] is the fundamental results that characterize the feasibility. 











Theorem 1. [17] Consider a unicast network setting G and assume I — G is irreducible. The 
following statements are equivalent: 

1) There exists a power vector p > 0 such that Gp > 0. 

2) A(I - G) < 1. 

3) G 1 = YZoi 1 ~ G) fc exists and is positive component-wise, with linifc_ s . 00 (I — G ) k = 0. 
Moreover, there exists p > 0 such that Gp = 0 if and only if A(I — G) = 1. 

For a multicast system, the main result of this paper is the following theorem. 

Theorem 2. Consider a multicast network setting A(p) and assume it is irreducible, i.e., the 
matrices I — G for G G Qip.) are all irreducible. There exists a power vector p > 0 such that 
A(/z)p > n(p) if and only if maxGe6(/x){A(I — G)} < 1. 

Theorem [2] basically says that for a wireless multcast system, an SINR vector p is feasible if 
and only if p is feasible to any of its embedded unicast system. 


Corollary 1. When there are only two multcast sessions, i.e., N 
determined by the unicast system specified by 

9 r k i t ■ 

where k* = arg max{/i,— L — L } 1 i = 
fciS/Ci g ki 

i ’ 4 

That is, p is feasible if and only if A (I — G*) < 1. 



= 2, the feasibility of p is 

1,2, j ^ i. 


Corollary |T| follows straightforwardly from Theorem [2] Note that when N = 2, for any G G Q, 

we have _ 

A (I — G) = , / pi —x P '2 ■ 


9 r k 2 
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(4) 


So max Ge6 {A(I - G)} = A(I - G*). 

The proof of the necessary condition for Theorem [2] is straightforward. Suppose there exists 
a power vector p > 0 such that Ap > n. Then for any G G G, we have Gp > n G > 0. By 
Theorem [TJ A(I — G) < 1 for all G, which implies max Ge g{A(I — G)} < 1. 

In the rest of this section, we prove the sufficient condition and assume that max Ge g{A(T — 
G)} < 1. By Theorem [T| it indicates that for each G G Q, G 1 > 0 exists, and thus p = 
G > 0 exists. For each receiver, define 

A = {pG M. n : a- p > nf*, p > 0 } . 








Note that af' is a row vector as defined before. A 1 '/ is an intersection of half-spaces and thus is 
convex. Our proof is based on Helly’s theorem given below. 


Theorem 3. (Helly’s theorem) [18]. Let T be a finite collection of com ex sets in M . The 
intersection of all the sets of T is non-empty if and only if any N + 1 of them has non-empty 
intersection. 


In our case, 

{Af = G/Q}. (5) 

There are in total K convex sets in T. If any N +1 of them have non-empty intersection, then all 
of them have non-empty intersection, i.e., the SINR is feasible. The number of all combinations 
of such iV +1 sets is ( V A , J. We first show the proof for N — 2. Then we use the mathematical 
induction to prove the general case. 

Lemma 2. Suppose X = (X i3 ) is an N x N matrix satisfying X l3 = 1 for i = j and X t] < 0 
for i f j. Let S be a subset of {1,..., N} and X' be the matrix by removing the i-th row and 
i-th column of X for all i G S. If A(I — X) < 1, then A(I — X') < 1. 

Proof: Since A(I — X) < 1, by Theorem [Tj there exists a vector p > 0 such that Xp > 0. 
Let p' G M ;V I 5 be the vector constructed by removing the i-th entry in p for all i G S. Since 

Xj.j < 0 for i f j , it can be verified that X'p' > 0, which implies A(I — X') <1. ■ 

Lemma 3. Consider G.G £ (? such that G differs from G only in one row. i.e., ki f ki for 

one i G {1,..., N} and kj = k 3 for j f i. Let p = G _1 riQ and p = G^riQ. There exists 

P £ {p, p} such that Gp > rig and Gp > n^,. 

Proof: Since p = G^ 1 !!^, and p = G _1 riQ, we have Gp = and Gp = n^. If p = p, 
automatically we have Gp = n c; and Gp = n^,. In the following discussion, we consider the 
case when p f p. Without loss of generality, assume that G and G differ in the first row, that 



is, k\ 7 ^ k] and k, = kj for j ^ 1 . Let us partition G into four blocks as follows. 
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Similarly, G is partitioned into four blocks as G = 


1 

B 

C 

D 


Note that G and G share the 

same three blocks: 1,(7 and D. We consider Gp — n G and Gp — n G . Since a."’ = a-'* and 
n h = n k, for a- p = afp = n k . = n*. and afp = afp = n*. = n*. for 

i = 2,.. ., N. If a^p = n ki or a^p = n kl , then Gp = n G or Gp = n G , which implies p = p. 
Therefore when p / p, we must have a, 1 p ^ n ki and a^ p ^ . In the following we prove 

that, either a^p > n k] or a^p > n ki but not both, that is, (a^p — n^Xa^p — n ki ) < 0. 
Since G _1 > 0 exists, by Theorem [l] and Lemma [l] I) 1 exists. By block-wise inversion 


1191, the inverse of G can be written as 

a — aAD ~ 1 
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Similarly, we have 


a ^"P ~n~ kl = -ba \BD 1 n'- n ki ) + {AD 1 n -n k J. 

Then 

(a^p-n^Xa^p-n^) 

= —a6 _1 (AZi)~ 1 n / — n~ ki ) 2 — ba~ 1 {BD^ 1 n — n~ ki ) 2 + 2(AD~ 1 n / — n ki ){BD~ 1 n' — n~ ki ) 

= — \Jab~ 1 {AD l n — n k] ) — Vba~ 1 {BD^ 1 n — n ki ) 

< 0 . 

Further, since a^p ^ n ki and a 1 p ^ n k , ( a 1 p — n kt )( a 1 p — n ki ) < 0. In summary, there 
exists p G {p, p} such that Gp > iIq and Gp > n<~,. ■ 


A. Two Multicast Sessions N = 2 

If K = 2, i.e., K[ = K -2 = 1, it is the unicast scenario and Theorem [2] is true straightforwardly. 
If K\ = 1, K‘> = 2 or I\\ = 2, Jl 2 = 1 or I\\ = K- 2 = 2, then any three subsets of J~ must be 
A) , A{. A kj for i = 1 or 2 and j ^ i. Let 
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a 1 
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and G = 


kn 

L a /J 

kn 

L a /J 


By Lemma 3} there exists p such that Gp > and Gp > n^, which implies p e {A} 


n 


A\ A A/). Further by Helly’s theorem, the intersection of all sets in T is non-empty. For other 
values of K\ and Ji 2 , we divide the { k1 + K '-) combinations of three sets of T into two parts: 
1) two sets belong to transmitter T, and one set belongs to transmitter Tj where j ^ L 2) three 
sets belong to the same transmitter Tj for i = 1 or 2. In the first case, the three sets could 


be A ki , A ]', Aj 1 . We use the same argument as before, and conclude that the three sets have a 

k y k" 

non-empty intersection. In the second case, the three sets could be A t \ A ^. It is easy to 

, y yt 

verified that p, = maxjrij } and pj = 0 is one of their intersection points. Overall, we 

prove that any three sets of T have a non-empty intersection, and thus the intersection of all 
sets is non-empty. 













B. Multicast Sessions with general N 


We use mathematical induction to prove Theorem |2j We already show that it is true when 
N = 2. Assume that the theorem holds for all numbers less than or equal to N — 1 and now we 
prove that it also holds for N. If K, = 1 for all i, it is the unicast scenario and Theorem [2] is 
true. Otherwise, we categorize the combinations of N + 1 sets of T into N parts: 1) Receivers 
of N transmitters are involved: A \', A k 2 2 , ■ ■ ■, A!1^' , A ^ 1 . 2) Receivers of A r — 1 transmitters are 
involved: A kl ,..., Afs^ , Af+i ,..., A 1 ^ 1 , A^, A^ where i,l ^ j. ■ ■ • D) Receivers of N — D +1 
transmitters are involved. ■ ■ ■ N) Receivers of 1 transmitter is involved. 

We prove the first part. Let 
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By Lemma 3 there exists p such that Gp > and Gp > n^, which implies p € (fl jLiA^ D 

We prove the D) part for D = 2,..., N. Suppose the I) — 1 transmitters whose receivers 
are not involved in the N + 1 sets, are di, d 2: ..., d D _i e {1, ■ ■ • ,N}. We simply let p dl = 
p d2 = • • • = Pd D _ 1 = 0. The resulting system is equivalent to having N — D + 1 multicast 
sections characterized by matrix A', which is constructed by removing the rows in A that 
corresponds to the receivers of transmitter d and the d -th column of A, for d = d\..... do- 1 - 
Define Q' C ]^( N - D + 1 ) X ( N ~ D + 1 ) for A'. For any G' £ Q', we can find a G £ Q such that, G' 
is constructed by removing the d-th row and d -th column of G for all d = d\,... ,dn- 1 - Since 
A(I—G) < 1 for all G £ Q, by Lemma[2j A(I—G 7 ) < 1, and therefore max G 'G6'{A(I—G')} < 1. 
By the inductive hypothesis, we can apply Theorem [2] when N — D + 1 < N, and thus there 
exists p' > 0 such that A'p' > n', where n' is obtained by removing the entries that correspond 
to the receivers of transmitter T d for d — d\,.... do -1 ■ By inserting 0 back into p' at the position 
of transmitter T d for all d — d \,..., do- i, we get a power p > 0 which is in the N + 1 subsets. 

Overall, we have proved that any A r +1 subsets of T has a non-empty intersection. By Helly’s 
theorem, all subsets in T has an intersection. This completes the proof of Theorem |2j 








IV. Feasible SINR region and algorithm 

In this section, we characterize the feasible SINR region of a wireless multicast system by 
analytically obtaining its boundary points. By Proposition [T] we know that the feasible SINR 
region is downward comprehensive. That is, if /x is feasible, then any /x' satisfying 0 < /x' < V 
is also feasible. Therefore, finding the boundary points is enough to figure out the feasible 
SINR region. Our approach is to find the farthest point from the origin in a given direction. In 
mathematics, the problem is formulated as 


sup (3 

p 

s.t. A(/3/x)p > n(/3/x) 

P > 0, 

where /x is a given direction. By Theorem [2j there is a feasible solution to the above problem 
if and only if 


max {A(I - G)} = (3 • max {A(I - G)} < 1. 

GG^(/3/x) GG^(/x) 


That is 
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max Ge6(M) {A(I - G)}' 


Therefore, the optimal value is 


PM = 


max Gee(/x) {A(I - G)}' 

/3*(/x)/x is a boundary point of the SINR region.The open line segment defined by {ctp : 0 < 
a < f3*(/j,)} is in the feasible SINR region T, but is not in the feasible region if a > /3 *(/j,). 

We note that the size of Q is fl^i Ku which grows exponentially with A. It is not an efficient 
method to calculate the Perron-Frobenius eigenvalue of all the embedded unicast systems and 
find out the maximum one. Next, we propose an iterative algorithm to compute /3*(/x). For 
i — 1, 2,..., N, let ej denote the iV-dimensional column vector such that the /-1h component of 
ej is 1 while the others are 0. The algorithm is described in Algorithm |T] 

For receiver /?,■’'*, (ej — &f l )p ( - k ' 1 is the sum of the interference power. The power of transmitter 
Tj is updated by the maximum interference power experienced by the receivers in its multicast 


session. This idea is similar to the distributed power control algorithm for unicast systems 201 


to solve the power balancing problem. Recall that in [20], given a normalized interference link 






Algorithm 1 Iterative algorithm 


l: Choose p (0 ) G W N > 0 and k 4— 0 

2: repeat 


3: 

4: 

5: 

6 : 

7: 


for i = 1 to N do 

Vi k) ma x k . eK . {(ef - af^pW} 

end for 

<— minf =1 { C 
(fc+l) , yW 

p ^ l|y (fe) l! 

k <— k + 1 


9: until convergence 

10 : return (3^ 


gain matrix I — G, the algorithm works as p( fc+1 ) = — j— where k is the iteration index. It 
is well known that when I — G is primitive (to be defined later), 11 (I — G)p^ 11 converges to the 
Perron-Frobenius eigenvalue of I — G, and p' fc) converges to the corresponding eigenvector. In 
our proposed algorithm, we are dealing with multicast systems. For notation simplicity, define 

Z(p) = {Z = I-G:Gg£(/x)}. (6) 


Z includes the normalized interference link gain matrices of all the embedded unicast systems 
and Z > 0 for all Z G Z. Given any p > 0, due to the structure of Q ( 71 ), there always exists 
Z G Z such that Zp > Zp for all Z G Z. Our algorithm works as p( fc+1 ) = , where 

Z W G Z is chosen such that > Zp^ for all Z G Z. In the rest of this section, we 

show the convergence of the algorithm. 


Lemma 4. The sequence {/3^} generated by Algorithm [7] is monotonically increasing and 


bounded above by 


^{A(I—G)} = max ze ^ ){ A(Z)} ’ and thus is convergent. 


Proof: By Algorithm [TJ we have y ^ > Zp^ for all Z G Z, and p^ > Then 


(*) 


z p'' M> = z A m < z 


P 


(fc) 


r( k ) 


< 


P 


(fc+1) 


/3( k ) | |y( fc ) 11 /)( fc )||y( fc )|| j3( k ) 














Since the above inequality holds for all Z G Z, we have 


/3 (fe) < min 

ZG2 


N 

min 

i=l 


P 


(fc+1) 


[Zp( fc+1 ) 



N 

min 

Z=1 


Pi 


(fc+i) 


.(fc+i) 


/ 3 (fe+1) . 


That is, is monotonically increasing. On the other hand, since > P^Zp^ 

for all Z e Z, that is (I — /3^ Z)p^ > 0, we have X(f3^Z) < 1 by Theorem 0 Therefore 
f j(k) < ma* zg ^ ){ A(z)} > and thus is convergent. ■ 

Denote lim^oo = [3*. Before we proceed to show that (3* = m . i% . z 1 { X(Z) } > we introduce 
the concept of primitive matrix and primitive set. 


Definition 4. [ 161 A square nonnegative matrix X is called primitive if there exists a positive 
integer n such that X” > 0. 


The class of primitive matrices is a subclass of irreducible matrices. If X is primitive, then 
its Perron-Frobenius eigenvalue is strictly greater than all other eigenvalues in absolute value. 
The primitive condition guarantees the convergence of the aforementioned distributed power 
control algorithm for unicast systems. In our multicast case, we need to use the concept of 
primitive set, which replaces a single matrix and powers of that matrix with a set of matrices 
and inhomogeneous products of matrices from the set. 


Definition 5. [211 Let Z be a set of N x N nonnegative matrices. For a positive integer n, 
let 0(n) be an arbitrary product of n matrices from Z, with any ordering and with repetitions 
permitted. Define Z to be a primitive set if there is a positive integer n such that every 0(n) is 
positive. 


It can be seen that a necessary condition for Z to be primitive is that Z is primitive for all 
Z G Z. One of the sufficient conditions for Z to be primitive is that for any Z e Z, in each row 


and each column of Z, there are more than half of the entries that are positive [21]. Interested 
readers can refer to [21] for more information of the primitive set. It needs to be mentioned that 
when the system is composed of two multicast sessions, Z e Z are always non-primitive, and 
therefore Z cannot be primitive. However, for this case, Corollary |T] already gives an explicit 
and simple solution to the feasible SINR region. Algorithm [T] works for systems with more than 
two multicast sessions and with Z being primitive. 










iteration step k 


Fig. 2: Convergence of Algorithm [j] In this example, there are four multicast sessions and each 
has three receivers. The link gains are randomly drawn from a uniform distribution on [0,1). 


Theorem 4 . If the matrix set Z(p ) defined in ([6]) is primitive, then (fi k) converges to 3* = 
maxz ^ ) { A(Z) | for an arbitrary initial value p (0) > 0. Moreover, p- ' converges to a power 
vector p* such that lim^t*, T(ap*) achieves the boundary point (3* p. 

The proof is provided in Appendix [A] Fig. [2] illustrates the typical behavior of Iterative 
Algorithm [T] In this example, there are four multicast sessions and each has three receivers. 
The link gains are randomly drawn from a uniform distribution on [0,1). As we can see, (3^ 
converges within a small number of iterations. By using Algorithm [T] we can efficiently check 
the feasibility of an SINR vector p by checking the value of (3*(p). If (3*(p ) < 1, p is infeasible, 
and vice versa. Besides, Algorithm [T] can be used to find the optimal solution of the classic power 
balancing problem for multicast systems, which is in the following form 

N 

sup min 7«(p) 

P i=1 

s.t. p > 0, 


and the solution is /3*(1). 


















V. Geometric Properties of the Feasible SINR region 


In this section, we discuss the geometric properties of the feasible SINR region. Let D(pf) 
denote the diagonal matrix constructed by 


£>(/*) 


!M 0 
0 jl2 


0 

0 


0 • • • 0 

By Theorem [2j the feasible SINR region is equivalent to 

T = Lg 


N : max (A(Z)} < 1 

Z 




D N 


max 

zgz(i) 


{a(D|»Z)}<i} 


= n {/* 6 R" : MD(m)Z) < l}. 

zez(i) 

That is, the feasible SINR region of a multicast system is the intersection of the feasible SINR 
regions of all its embedded unicast systems. Let Y c = \Y denote the complement of T in K+, 
i.e., the infeasible SINR region. Next, we investigate the convexity of Y c and the log-convexity 
of T. 


A. Convexity of Y c 

For unicast systems, it has been proved in [5J that the infeasible SINR regions of a general two 
user system and a general three user system are convex. It is also shown in |4] that the convexity 
of the infeasible SINR region does not hold for a general four user system. For multicast systems, 
we have the following observation. 

Theorem 5. The infeasible SINR region of a general system consisting of two multicast sessions 
is convex. The convexity property does not hold for a general system consisting of more than 
two multicast sessions. 


When there are two multicast sessions, by Corollary [I] the feasible SINR region is 


Y — < [/ii,/i2] G ®Fj_ : pip2 < 


Ok * Ok* 

rY.fi *r 2 2 ,t 2 

Ok* Q k% 

rf,t2 r 2 2 ,ti 






It is ready to verify that Y c is convex. 

When there are more than two multicast sessions, Y c is the union of the infeasible SINR 
regions of all the embedded unicast systems and is in general non-convex. Fig. [3] illustrates the 
Y c for a system consisting of three multicast sessions, where the link gain matrix is given by 
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It can be seen that its Y c is non-convex. 

B. Log-convexity of Y 

We first introduce the notion of log-convexity. Let log(//) = [log/ii, log/x 2 , • • •, log/ijv] and 
log(Y) = (log(p) : n e Y}. We say a set Y is log-convex if log(Y) is convex. Since log(-) : 
Y —y log(Y) is a bijective mapping, we have 

log(T)= p| { log(/t) e : A(£>(/t)Z) < 1 j. 
ze.z(i) 

It has been proved in |3} that the feasible SINR region of a unicast system is log-convex. So 
log(Y), the intersection of the SINR regions of all its embedded unicast, is also log-convex. We 
conclude this by the following theorem. 

Theorem 6. The feasible SINR region of a multicast system is log-convex. In other words, the 
feasible SINR, expressed in decibels, is a convex set. 

VI. Feasibility of SINR with Power Constraints 

So far, we have discussed the feasibility of SINR for a multicast system in the case of unlimited 
power. In this section, we consider that besides p > 0, the power vector are also subject to the 
linear constraints 

Y Pi 1,..., M, 





Fig. 3: Infeasible SINR region of a three-multicast-session system. 

where Q m C {1,..., N} and M is the number of constraints. When Q m = {1,..., N}, it is 
a constraint on the total power. When Q m = {z}, it is a constraint on the individual power of 
transmitter Tj. Define the power set by 

P = {p>Oand = 

Now the feasibility of SINR vector p is decided by whether there exists p G V such that 
A(p)p = n(p). Note that the power vectors in V are downward comprehensive. That is, if 
p' G V, then p e "P if 0 < p < p'. Hence using the same argument as in Proposition [I] we 
know that p is feasible if and only if there exits p £ V such that A(p)p > n(p). Our results 
generalize the feasibility condition derived in [[6] for a unicast system to a multicast system. 

Definition 6. For a matrix X G W KxN , a vector y G and a set VL C {1,..., N}, ^(X. y, 0) 
is the operation to add y to the j-th column of X, for all j e Q. That is, Z = p;(X. y, Q), where 
Zij = Xij + yi for all i G {1,..., K} and j G 0, and Zij = X, :l for the else. 

Theorem 7. Consider a multicast network setting A(p) and assume the matrices I — G for 
G G G(p) are all irreducible. There exists a power vector p G V such that A(p)p > nfp) if 




and only if 


max max 

Gee(/x) 


{ A (v,(I- G ,|=,Q m ))}< 1 . 

V Jrzlm ) 


Proof: It is already known from [J6] that, for a unicast system G, there exists p e V 
such that Gp > n G if and only if max me { lv ..m{A('0(I -G,£?-,n m ))} < 1. We first prove 

5 ’ P^m 

the necessary condition. Suppose there exists p e V such that A(/x)p > n(/x). Then for any 
G e Q(n), Gp > n G , which implies max^n M} {A(0(I-G, < 1. Regarding all 

yitm 

G e G(n), we have max G60(fl) max^n m{A(^(I - G, Q m ))}leql. 

Film 

Next we prove the sufficient condition. For any G 6 G ([*-), since 0 < I — G < -0(1 — 


Gfor all m, by the Perron-Frobenius Theorem for irreducible matrices [ 16], A(I — 


G) < A(0(I — G, A 12 -, f2 m )) < 1. This implies that G 1 exists and p = G 1 n G e V. The rest 

PGlm 

of the proof follows the same argument as the proof of Theorem |2j ■ 

Note that 


max max 
GeG(n) me{i,...,M} 


n G 




{ A(Wl — G, 1 = max max |a(V(I — G, 

l V v p Qrn V J me{i,...,M}Ged(j*)\ V v 

Similar to (|6]), for each of the M linear constraints, define 

2t.M = U-GA,(l„):G€5M}. 

I Pn m J 

By using Algorithm [I] with Z< lin (ji), we can find a supremum /3^ m (n). The farthest point of 
the SINR region in direction /i is then min^ =1 {pQ rn (p) }p. By this approach, the feasible SINR 
region is characterized. On the other hand, if > 1, // is feasible. Fig. |4j plots 

the feasible SINR region of the network example in Fig. [T] with a power constraint on the total 
power. In this example, the link gain matrix is 

R\ 

K\ 

R\ 

Rl 

and the power constraint is p\ + p 2 < 2. The four dashed lines are the boundary of the feasible 
SINR regions of four embedded unicast systems and the solid line is the boundary of the multicast 
system. It can be seen that under power contraint, the infeasible SINR region is not necessary 
to be convex even for a multicast system with two multicast sessions. 
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Fig. 4: Feasible SINR region for a multcast system with two multicast sessions, under a total 
power constraint. The dashed lines correspond to the four embedded unicast systems and the 
solid line corresponds to the multicast system. 


In the end of this section, we introduce an application of our multicast model to a time varying 
unicast system. Consider a unicast system consisting of N transmitter-receiver pairs, where the 
channel gains among them vary with time due to the mobility of the receivers. Let h*(f) for 
i = 1,... ,N denote the link gain vector from N transmitters to the z-th receiver at time t. As 
argued in [6], h ( (7) can be modeled with discrete states, that is, h,;(t) is randomly selected from 
a Ante set {hf, hf,..., hf ; } for all i. An SINR vector fx is said to be zero-outage if there exists 
a power such that no matter what the link gain realization is, the SINR is achievable all the time. 
Such a zero-outage SINR problem can be mapped to a feasible SINR problem of a multicast 
system. The idea is to let one receiver II, pretend to be K, receivers, i.e., R},..., Rf\ and /?,f' 
only experiences the link gain h . The feasible SINR region of this artificial multicast system 
is exactly the zero-outage SINR region of the original time-varying unicast sytem. Theorem [2] 
and Theorem [7] can be applied. It needs to be mentioned that the scenario considered here is 
different from that in [16]. In [|6], the power can change with the channel states and is subject to 
an average power constraint. In our model, the power is universal for all possible channel states. 
















VII. Conclusion 


In this paper, we characterize the feasibility condition of an SINR vector for a multicast 
system, which generalizes the Perron-Frobenius Theorem for a unicast system. We also propose 
an iterative algorithm which can efficiently check the condition and compute the boundary 
points of the feasible SINR region. According to the earlier mentioned Gaussian interference 
assumption, by describing the feasible SINR region, we directly obtain the feasible rate region 
by applying the Shannon Capacity formula for AWGN channels that maps the SINR to the rate. 


Appendix A 
Proof of Theorem [4] 

Proof: Let Z (k) denote one of the matrices at the /, -th iteration such that Z (k> p ik) > ZpW 
for all Z £ Z. From the construction of the algorithm, we have 


Z ( V fc) < -4rrP (fc) for all fceN. 


(7) 


Moreover, there exists 1 < i < N such that [Z^p( fc )]j = 7 ^- [p' 7 '"*],;. We note that each vector 
p (fc ) is a unit vector, as |p7 ,:) | = 1. By the Bolzano-Weierstrass Theorem and the compactness 
of the unit ball in there exists a convergent subsequence, that is, p !/i: L —>• p*. By Lemma 
|4j — y P*. Suppose at p*, Z* e Z is one of the matrices that satisfy Z*p* > Zp* for 

all Z e Z. Taking the limit of ([7]) with respect to the subsequence indexed by k v we have 
Z*p* < 4rP*. If Z*p* = tLp*, since Z* is irreducible, /3* = > 


other hand, (3* < 


A(Z*) — max zezM {\(Z)} 

by Lemma 4 Therefore, P* = 1 


. On the 


max zg ^ ((i ){A(Z)} J |_f ’ M max zez{(i) {A(Z)} ‘ 

If Z*p* f -j-p*, since Z is primitive, there exists integer n such that an arbitrary product 
of n matrices from Z is positive, i.e., 0(n) > 0, and therefore 0(n)Z*p * < 0(n)^p*. By 
the continuity of the mapping, there exists p^ close enough to p* such that 0(n)Z*p( fc ) < 
0(/■;,) 7-p i:/ '^ and Z' :k) = Z*. We now apply the algorithm for n more iterations from \f k) . For 
i = 0,l,...,n — lwe have the following inequalities 

2(fc+*+i)p(fe+0 < 2( fc +*)p( fc +*) (8) 

due to that the selection matrix satisfies Z^ k+ ' i:> \f k+l ) > Zp( fc+ 7 f or a \\ z e Z. Meanwhile by 
Algorithm [I] 

(k+i) _ 


P 


y(fc+i-l) r^(k+i-l)-p{k+i-l) Z (k+i- 1 ) . . . Z( fc+1 )Z^7p( fc ) 


T (k+i- 1) I 


|2(fc+*-l)p(fe+i-l) 


|Z(fc-H—i) . . . Z( fc+1 )Z( fc )p(0 












By substituting p( fc+ *) into ([8]), we get 


z (k+i+l) z (k+i— l)r^(k+i— 2) _ _ _ 20+l) Z 0)p(k) < Z {k+i) z (k+i— l)r^(k+i— 2) _ _ _ 20+l) Z (k)p(k) (9) 

Let us take a look at these inequalities step by step. By ([8]) for i = 0, Z^' +1 )p( fc ) < Z^'p^. By 
multiplying Z^ k+ ^ on both side of the inequality, we have 

2(fc+2)2(M-i)p(fc) < 2^ +2 ^Z^^p^ (10) 

By ([9]) for i = 1, Z^'+^Z^p^ < Z( fc+1) Z( fc )p( fc ). Along with ( [TO] ), we have 

2( fe+2 )z( fe+1 )p( fc ) < Z^+^Z^p^. 

By multiplying Z^ fc+3 ^ on both side of the above inequality, we have Z^+^Z^+^Z^+^pO^ < 

Z (fc+3) z (fe+i)2( fc )p( fe ). By @ for % = 2, z^+^Z^Z^p^ < z^ k +^Z^ k+1) Z^p ik l So 

2(fc+3) z ( fc + 2 )2( fc+1 )p( fc ) < z ( fc+2 )z^ +1 -*Z^^p^ 


By repeating this procedure for n — 1 times, we can finally get 

Z ( fc+n ) z ( fc+n_1 )... Z ( fc+1 )p( fc ) < z l ' k+n - 1 )Z ( ' k+n - 2 ') ■ ■ ■ Z^^Z^p^ 


(ID 


Since 0(n)Z*p® < 0(n)^rp^ holds for arbitrary @(n), we let @(n) = z < ' k+n ' > Z ( ' k+n 0 ■ ■ • Z^ k+1 \ 
Along with ( fTT| ), we have 


z (k+n) z (k+n- 1) . . . Z ( fc +0Z*p (fc) = 0(^)Z*p (fc) 

< e(n)kp<‘) 

1 


— _ z {k+n) z (k+n-l) _ _ . Z (k+1) (k) 

P* ^ 


By multiplying 


| 2 (fe+n-l)... 2 (fc+l)z*p( fe )| 


Z (k+n)^(k+n) _ z (k+n) 


_ z (k+n—l) z (k+n—2) _ _ _ Z^^^ Z* 

“ P* ” ’ P 

on both side of the inequality, we have 

Z (k+n-l) . . . Z ( fc + 1 )Z*p( fc ) 


Z (fc+n—1) . . . Z ( fc + 1 )Z*p0)| 

1 iPk+n-l) . . . Z( fc + 1 )Z*p^ 





This implies / 3 ( ' k+n > > ( 3 *, which contradicts with that ( 3 * is the limit. Hence there must have 
P — ^rP • 

We prove p ^ —>■ p* by contradiction. Suppose there exists another subsequence such that 
p k i —y p' and p* ^ p'. Then Z*p' < ip'. Meanwhile we already have Z*p* = 4rP*. By the 


Subinvariance Theorem in [16] (pp. 23), p* = p\ which contradicts with the assumption that 
p* 7^ p'. Therefore p (/ ^ converges to p*. 

Since Z*p* = -^p* and Zp* < -^p* for all Z <G Z{n), it is ready to see that lim^oo 7j(ap*) = 
min Ze2 (i){p^|-} = for all i = 1,..., N. So lim^^ T(ap*) = /T>. ■ 
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